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Abstract 

We consider large- multi-matrix models whose action closely mimics that of 
Yang-Mills theory, including gauge- fixing and ghost terms. We show that the fac- 
torized Schwinger-Dyson loop equations, expressed in terms of the generating series 
of gluon and ghost correlations G(^), are quadratic equations S^G = GS^^G in con- 
catenation of correlations. The Schwinger-Dyson operator 5* is built from the left 
annihilation operator, which does not satisfy the Leibnitz rule with respect to con- 
catenation. So the loop equations are not differential equations. We show that left 
annihilation is a derivation of the graded shuffle product of gluon and ghost corre- 
lations. The shuffle product is the point-wise product of Wilson loops, expressed 
in terms of correlations. So in the limit where concatenation is approximated by 
shuffle products, the loop equations become differential equations. Remarkably, the 
Schwinger-Dyson operator as a whole is also a derivation of the graded shuffle prod- 
uct. This allows us to turn the loop equations into linear equations for the shuffle 
reciprocal, which might serve as a starting point for an approximation scheme. 
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1 Introduction 

QCD, a quantum Yang-Mills theory is the best candidate for a theory of strongly interacting 
subatomic particles. It is an outstanding challenge to understand its non-perturbative features 
and their mathematical formulation. Yang-Mills theory includes as dynamical degrees of free- 
dom, a collection of N x N matrices (gluons), where = 3 in nature. The limit as the number 
of colours ^ oo is a promising starting point to study this theory [1]. Large- matrix models 
[2] are simplified models for the dynamics of gluons. 

In this paper, we establish some properties of a class of large- multi-matrix models that 
may be regarded as toy-models for gauge-fixed Yang-Mills theory. These models have both her- 
mitian complex matrices (bosonic gluons) and hermitian grassmann matrices (fermionic ghosts) 
and we call them Yang-Mills matrix models with ghosts. However, these models are generically 
not supersymmetric. Our work is inspired by that of Makeenko and Migdal jSj |4l [5] on the 
loop equations. Matrix models and their loop equations may be formulated in several different 
ways [U El El [HI [ini [m [I2j. There is a large literature on matrix models involving both bosonic 
and fermionic degrees of freedom, especially in the context of the matrix approach to M-theory 
|13l [T4l [TS] and in the matrix regularisation of the supermembrane |16j . 

'Solving' a matrix-model can be regarded as determining the gluon and ghost correlation 
functions. We obtain quantum corrected equations of motion (large- A^ factorised Schwinger- 
Dyson or loop equations) for these correlations. They involve a 'classical' term linear in correla- 
tions, coming from the variation of the action. This is the 'Schwinger-Dyson operator' 5* acting 
on correlations. 5* is built from the left annihilation operator. The loop equations also involve 
a quadratic 'quantum' term in correlations, coming from the change in path integral measure. 
This involves concatenation products of gluon and ghost correlation tensors. However, the left 
annihilation operator does not satisfy the Leibnitz rule (i.e. is not a derivation) with respect 
to concatenation. So the loop equations are not differential equations. This is perhaps part of 
the reason why the loop equations have been difficult to solve, though they were derived for 
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the Wilson loops of Yang Mills theory over 25 years ago. It is therefore imperative to uncover 
any hidden mathematical structures of the loop equations, which may help in solving them and 
placing them in their natural mathematical context. 

The great success of calculus in solving problems of classical mechanics is due to the fact 
that the equations are differential equations, rather than, say, difference equations. So it is 
interesting to know whether there is some limit or approximation where the loop equations 
become differential equations. The main result of this paper is that this is indeed the case, 
in the limit where concatenation of correlations is replaced by their shuffle products. While 
concatenation arises from concatenation of loops, shuffle arises from the point-wise product 
of Wilson loops. We show that this picture is robust, in the sense that it is not spoiled by 
the inclusion of gauge fixing and ghost terms in the action. More precisely, we show that left 
annihilation as well as the Schwinger-Dyson operator 5* , satisfy the Leibnitz rule with respect 
to the (graded) shuffle product of gluon and ghost correlations. The latter allows us to reduce 
the non-linear loop equations to linear equations for the shuffle-reciprocal of the generating series 
for correlations. Though this is not a property shared by generic matrix models, it does carry 
over to 3 + Id gauge-fixed Yang-Mills theory. For, all we use is the algebraic structure of the 
action, and general properties of the large- limit. 

In a previous paper [12j, we proposed an approximation scheme to compute correlation 
functions by solving the loop equations, in the context of bosonic matrix models. We expanded 
the concatenation product around the shuffle product and used the shuffle reciprocal to reduce 
the loop equations to linear differential equations in the shuffle algebra at the 0**^ order of 
the expansion. In simple cases, this was shown to give a rough approximation to the exact 
correlations. Aside from this practical application, a mathematical lesson from our work is that 
it may be fruitful to view the loop equations as living in the differential bi-algebra formed by 
shuffle, concatenation and their derivations. 

The property that 5* is a derivation of the shuffle product, is a finite and differential- 
algebraic reformulation in terms of correlation tensors, of a property of large- Yang-Mills 
theory mentioned in [4j. 5* is analogous to the path derivative of the area derivative operator, 
which was said to satisfy the Leibnitz rule with respect to point-wise products of Wilson loops. 
We hope that our alternative viewpoint and finite formulation of these infinite dimensional 
notions is useful to better understand the loop equations. 

2 Gluon-Ghost correlations in the large- limit 

Motivated by the Lagrangian of gauge-fixed Yang-Mills theory (say in a class of covariant gauges 
labeled by 

C = tr |i9^^,(9^A'^ - d'A^') - igd^,A,[A^ , A^] - ^[^^, 

+ Y^{d^A^f + d^cd^c - igd^c[A^,c]Y (1) 

we consider models with A matrices Ai, 1 < i < A, which are N x N matrices in 'colour' 
space. We will call them Yang-Mills matrix models with ghosts if their action is of the form 

tr S = ^ tr C'^AiAj + tr C'^''Ai[A,,Ak] - ^ tr [A, A,][Ak, Aijg^'g^'. (2) 
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We call all matrices Ai irrespective of whether they are gluons {Af^), ghosts (c) or anti-ghosts 
(c). Indices i,j,k are shorthand for position coordinates and polarisation indices. They also 
specify whether the matrix is a gluon, a ghost or an anti-ghost via the ghost number of an index 



0, if A j is a gluon; 

1, if Ai, is a ghost J- . (3) 
— 1 if j4j , is an anti-ghost 



If t° are hermitian generators of the Lie algebra of U{N) , then Ai = Aft^ with Af* = Af . 
For a gluon, ^{i) = and Af is a real number while for a ghost or anti-ghost, = ±1 and 
Af is a real grassmann variable. Moreover, 

{AfAfy = (4) 

Note that the ghost matrices are not related to the anti-ghost matrices by hermiticity. 

It is also useful to define the ghost number of a tensor (7*1 ■ to be that of the multi- index: 
#(^1 ■■■in) = Z]fc=i#(^fc) 0- In keeping with the structure of ([T|) we allow gluon, ghost and 
anti-ghost matrices in the quadratic and cubic terms of ([2]) but only equal numbers of ghosts 
and anti-ghosts in any term. In the cubic term, [ylj, ^fc] denotes the anti-commutator if neither 
Aj nor ylfc is a gluon and the commutator otherwise. In the quartic term of ^ we allow 
only gluons. In other words, we assume g^^ vanishes if either i or j corresponds to a ghost or 
anti-ghost index and that the ghost numbers of C*-' and C^^^ vanish. ([1]) can be regarded as a 
limiting case of ([2]) for appropriate integral kernels C^^,C^^^ and g^^ , when the indices become 
continuous. 

The action of our matrix model will be written as tr S{A) = tr S^Aj . It defines coupling 
tensors . The partition function is defined as Z = J IljdAje~^ ^^^^ where the integration 
is over all independent matrix elements of the gluon, ghost and anti-ghost matrices. Observables 
are correlation tensors 

C^A, . . . = i / n,.A,e- " -<-)|^, ^ ^ ^ '^A,,. (5) 

They are symmetric in the 's up to a possible sign. For example 

l'^A,'^A,) = (-lf<->*'^H'^A,'^A,). (0) 

This is because matrix elements [^i]^ are graded commutative. So we pick up a minus sign 
under transposition of ghost or anti-ghost matrices 

[^i]?[^.]S = (-l)*^^)*^^'^[^i]S[^.]^ (7) 

It follows that the trace operation is graded cyclic, for example tr AiAjA^ = (— tr A^AiAj . 

We restrict to actions Aj whose non-vanishing coupling tensors have zero ghost number. 
Gauge-fixed Yang Mills theory ^ is of this type. In such a theory, correlations of a tensor with 
non-zero total ghost number must vanishU 

i^Aj, . . . ^^Aj„) = if + #(/2) + • ■ • + #(/„) / 0. (8) 



^For a general tensor G = Gi(,^ , i^{G) = n if #(/) = n for each I for which Gi is non-vanishing. A 
general tensor has a well-defined ghost number only if all terms have the same ghost number. Capitals denote 
multi-indices I = iii2 ■ ■ ■ in and repeated indices are summed. 

^Consequence of Feynman rules: interaction vertices have equal numbers of ghost and anti-ghost lines attached. 
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For example tr AgAgA^Ac) = but tr AgA^Ac) can be non-trivial, where g,c and c 
stand for gluon, ghost and anti-ghost. Multi-trace correlators factorise into single trace correla- 
tions in the large- A*" limit: 

Gi= lim i^Ai), e.nd {^Aj,-..^^Aj J = Gj, ■ ■ ■ Gi„ + 0{1/N'). (9) 

Ai— >oo iV iV iV 

The hermiticity (jj]) of the matrices in colour space implies an order reversal property under com- 
plex conjugation G*^^^...^^ = {-iyGi„i^_^...i^i^ where p = #(ii)#(i2 • • • i„) + #{i2)#{h ■■■in) + 
• • • + #(«n-i)#(^n) • The factorised correlation tensors Gj are in general graded cyclic: 

G,. = Gu = {-l)WW)Gji. (10) 

Similarly, the only part of that contributes is its graded cyclic projection 

S^' ^ -^[S^' + 5''^(-l)#(^)#W + •••]. (11) 
bo we assume that S^'^ = 5-^^(-l)#(^)#(-^) for all I, J, i.e. that is graded cyclic. 



3 Schwinger-Dyson to loop equations in the large- limit 



The loop equations are quantum corrected equations of motion for correlation tensors. To derive 
them, we consider changes of variable (vector fields) 



Ai^ A'i = Ai + vUi. 



(12) 



Here v( could either be a small real number or a grassmann number with ghost number #{vj ) = 
±1. However, we cannot change a complex number by a grassmann- valued quantity and vice 
versa, so we require that v( are such that 

#« = #K') + #W- (13) 
An example of such a change of variable is a BRST transformation in Yang-Mills theory 



{d>^Al)\. 



(14) 



Here vl cx A is a constant grassmann quantity with ghost number —1, i.e. #(i'/) = —1 and it is 
easily seen that the conditions #(i) = #{vl) + are satisfied for the BRST transformation. 

Under ([H]), Ak ^ Ak + 5^^'^vIAlim- So if we define the 'loop' variable ^i^- = tr , 



g-Af tr S-'Aj 

'dA' 



n 

p=i 



-N tr S-'Aj 



det 



dAi 



det 



d[A,]t + vi[Ai]t 



d[A 
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= l + iVV'5}''^'(-l)*^*^*^^'^^/i^/2 + •••■ (15) 

The sign in the Jacobian comes from moving the grassmann (left) derivative through the product. 
Requiring the invariance of (0) under the changes of integration variables p2|) leads to the finite- 
Schwinger-Dyson equations (SDE) 

viS-^^^'H<^j,jj,) = z;/^^^^(-l)#W#(^0($,^$,^) + _L ^ di^;''-v{{^L,jM,). (16) 

p=i 

In the large- limit, the factorised SDE ignoring the last term are 

vjS'^^'^Gj.ij, = ^^^(-l)#(^)#(^i)G,,G,,. (17) 
Since correlations G/^ vanish if #(/i) 7^ we get 

yigJiiJ2Qj^^j^^yhihQ^^Q^^ for all vector fields v. (18) 

Now taking v( to be non-vanishing only for a fixed i, I gives us the loop equations (LE) 

S^''^'Gj,ij, = 6'/'''Gi,Gi, (19) 
Using graded cyclicity of and G/ we write this as (|/| is the length of the multi-index /) 

\Ji\S'^'Gji = Sj^'^^Gi.Gi^. (20) 
The loop equations have the same form as for a bosonic matrix model (see |12]). 

4 Loop equations in terms of concatenation and left annihilation 

Define the left annihilation operator acting on the generating series of correlations G{^) = Gi^^ 

[D,G]i = G,i and • • • D,,G]i = [D,„...,,G], = G,,..,„,. (21) 

The ghost number ^(DjG) = #{j) if G(^) has zero ghost number. In terms of the concatenation 
product 

[FG]i = dj'^FjGK; F{i)G{i) = FiGj^', (22) 
the LE take the same form as in the purely bosonic case |12j : 

J2in + l)S^'-'"'D,„...,M0 = GiOeCiC), or S'GiO = GiOmO 

n>0 

where 5* = ^(n + l)SJ'i -^"*L>j„...j,. (23) 

n>0 

Both LHS and RHS have ghost number #(i). 

Graded commutator To simplify notation below, we introduce the graded commutator 

[D,,D,] = D,D, - (-l)#«#(^')i?,A. (24) 

It is graded symmetric [Di,Dj] = — (— l)**^*)*^-'^ [Dj, Dj] , and reduces to the commutator if 
either i or j is a gluon and to the anti-commutator if neither is a gluon. We use the same 
notation for commutators and anti-commutators of the Ai . 
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5 Example: Ghost terms in Yang-Mills action 



As an example, consider a class of matrix models inspired by the ghost terms in the gauge-fixed 
Yang- Mills action([T]) d^c d^c — ig d^c [A^^, c] . The action is 

tr 5 = ^ tr &^AiAj + tr &^^Ai[Aj,Ak] (25) 



where [, ] is the graded commutator (j24p . If all the matrices are gluons, then this is a zero 
momentum Gaussian + Chern-Simons-type matrix model. Here we allow gluon, ghost and 
anti-ghost matrices so that for special choices of C*-' and C^^^ , this action can also model the 
terms d^c d^c and igd^c [A'^,c] appearing in the gauge-fixed Yang-Mills action. We obtain the 
coupling tensors S^^ and S^^^ , which can be chosen to be graded-cyclic. Then we obtain the 
factorised LE and the differential operator 5* . First, write the action as 



tr 5 = ^ tr C'^AA, + tr &^>'A{AjAk - (-1)*^^^*^''^ A^A,} 
= i tr C'^A^Aj + tr Ic'^'A^k - {-lf^'^*^''^C"'^ 



(26) 



Thus the coupling tensors are 



Sij = ic^ and S'^'' = C'^'' - C'''^ . (27) 

From graded cyclicity of the trace, it follows that the coupling tensors are graded cyclic 

Qij ^ ^_-^^m#U)sji and S'^'' = S'^'K (28) 

As for the differential operator 5* , from (j23p we have 

= C^Wj + SS^'^'Dkj = C^'Dj + 3(^0^''' - {-l)*^^'^*^^a^''^DkDj. (29) 

The interesting question (which we will answer in the affirmative in section [8]) is whether 5* is 
a derivation of the graded shuffle product of ghost number zero tensors. Mere graded cyclicity 
of S"^^^ is not sufficient for this. Rather, it is useful to write 5* in terms of graded commutators 
of left annihilations. The reason is that the left annihilation Di is like a first order differential 
operator (vector field) when acting on the shuffle algebra. While products of vector fields DiDj 
are no longer vector fields, their commutators [Di,Dj\ continue to be vector fields and therefore 
define derivations of the shuffle algebra. 

In terms of the graded commutator (j24p . the action may be written as 

tr S = tr ]^C'^AiAj + tr &^^Ai[Aj,Ak] = tr ]^C'^ AiAj + tr &^''[Ai, Aj]Ak. (30) 
Graded symmetry of the graded commutator in turn implies graded symmetry of C"^^^ 

More precisely, even if these quantities were non-vanishing, they would not contribute to the 
action. Therefore they can be taken to vanish without loss of generality. Using this property 
twice we can write 
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= &'Dj + ?,&^'DkDj - ?,{-l)*^^'^*^^^a^'DjDk = a'Dj + ?,&^^[Dk,Dj]. (32) 

Thus 5* is a linear combination of the left annihilation and its (anti-) commutators. Finally, the 
factorised LE for the action (j25p are 

a'DjGii) + 3C'^^[Dk, Dj]Gi() = Cm'GiO- (33) 

More generally, 5* for the full gauge-fixed Yang-Mills action is a linear sum of the following 
combinations of left annihilation operators 

Dg, D,, D,, [D,,D,U , [D,,Dg]., [D,,Dg]., [Dg„Dg,]., [Dg„[Dg„Dg,U.. (34) 

Each of these combinations arises as the variation of one or more terms in the gauge-fixed 
Yang- Mills action ([1]). Dg comes from varying terms involving two derivatives and two gluon 
fields (e.g. ^{df'Af.y and ^f,A„{^''A^' - ^^'A'')). and come from df,c d^'c. The anti- 
commutator [Dc,Dc]+ arises from the variation of gd^c [A^,c] with respect to the gluon. The 
commutators [Dc, Dg] and [Dc, Dg] arise from varying the same term with respect to an anti- 
ghost or a ghost (as in the example above). [Dg-^,Dg^] originates from varying the term linear 
in derivatives: gdfj,A^[A^, A'^]. Finally [-D^i, [Dg^,Dg^]] has its origin in the term independent 
of momentum ^ [^^, A,,] [A^^, A""] . 

We showed in that 5* for the purely gluonic part of the Yang-Mills -|- Chern-Simons -|- 
Gaussian action is a derivation of the shuffle product of gluon correlations. Here we show that 
even when gauge fixing and ghost terms are included, 5* is a derivation of the graded shuffle 
product of gluon-ghost correlations. 



6 Graded shuffle product 

We will call the extension of the shuffle product (see [12\ [T7] ) to correlation tensors of gluon and 
ghost matrices by the name graded shuffle product. It is essentially the shuffle product with 
minus signs when ghost or anti-ghost indices are transposed. The definition is 

[FoG]i^ J2 {-ir^'''''''^Fi,Gi,. (35) 

hUl2=I 

Ii\Jl2 = I is the condition that Ii and I2 are complementary order preserving sub-strings of I. 
In other words, we riffle-shuffle the card packs Ji and l2- For example if / = 21^2^3^ then one 
permissible choice is Ii = and I2 = iiis while Ii = 1^12-, I2 = HH is not allowed. We call 
7(1; /i, I2) the ghost crossing number of the ordered triple (/; Ii, I2) . It is just zero for bosonic 
matrix models. More generally, the string I is transformed into the string I1I2 by a minimum 
number of transpositions of neighboring indices. Each transposition ipiq ^ iqip contributes 
#{ip)i^{iq) ■ #(v)'?T^(*<?) if both ip,iq are gluons, if both are ghosts or anti-ghosts and 
— 1 if one was a ghost and the other an anti-ghost. The sum of these contributions is the ghost 
crossing number. For example let 

I = iii2hiii^; h = h = 12^- (36) 

The sequence of transpositions may be 
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^ i,i,i2i5i3(_l)#{M)#fe)+#{M)#fe)+#fe)#fe) 

^ ili4^5«2i3(-l)*^'*''^*^''''^- (37) 

The gf/iosi crossing sign (— l)'i'(^;^i.^2) independent of the choice of sequence of transpositions, 
so we don't need to stick to the minimum number to find the sign. Moreover, 

(_l)7(/;^,i^) = (_i)7(/;^f,^)+#(i)#W. (38) 
In the sequel we wih be interested in tensors Gj of zero ghost number, since the others vanish. 

Preservation of zero ghost number: The graded shuffle product of two tensors of zero ghost 
number is again a tensor of zero ghost number. Suppose Ff^ and Gj^ have zero ghost numbers, 
= #(/2) = 0. From ^ and ©, if / = /i U Is , then #(/) = #{h) + ^{h) = 0. So 
(F o G)j has zero ghost number. 

Commutativity of graded shuffle: The graded shuffle product FoG is commutative if either 
F or G has zero ghost number: 

[FoG]i = J2 {-iy^'''''''^Fj,Gi,= {-iy^''''^''^+*^''^*^''^Fj,Gi, 

hUl2=I IlUh=I 

hUli=I hUl2=I 



7 Left annihilation is a derivation of graded shuffle product 



Di is a derivation of the graded shuffle product of two ghost number zero tensors: Di{F oG) = 
{D,F) oG + Fo (DiG) if = = 0. To show this, we write 

mF o G)]i = [F o G]u = Y {-lV^'''''^''^Fj,Gj,. (40) 

hUl2=iI 

Now either z E /i or i E /2 1 so 



mF o G)]j = {-lV^'''''''''^Fu,Gi, + {-iy^''-^'''''^^Fi,Gu2 

hUl2=I - 

= E \{-IV'^"'"''''\d,F]i,Gi, + (-1)^(^^'^^'*^^)f,JAG],, 



(41) 



Since i doesn't cross any index, (-l)7{»^;i^i,^2) = (_i)7{^;^i,/2) , Similarly (-l)7(»-f;-fi>*-f2) 
(_l)#W#{^i)+7a;/i,^2)_ Thus 



[A(F o G)]i = Y {-ir^'''''''HF>iF]i,Gi, + (-l)#«#(^i)+^(^;^i'^2)F,jAG]7, 

IlUl2=I 



(42) 



Since correlation tensors vanish for non-zero ghost number we can take = 0. Hence 

[Di{F o G)]i = [{D,F) o G]i + [Fo {DiG)]i. Thus A is a derivation of the shuffle product of 
two tensors provided each has zero ghost number. More generally, if no assumption is made on 
the ghost number of F and G, then 

[D,{FoG)]i = [{D,F)oG]i+ Y {-l)*'^^*^''^^^^'''''''^Fi,[D,G]i,. 

hUl2=I 



(43) 
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Graded commutator of left annihilation: In the bosonic theory, a commutator [Di,Dj] 
of derivations is a derivation of the shuffle product [12]. More generally, we will show that if 
F, G have zero ghost number, then the graded commutator of left annihilations is a derivation 
of their shuffle product: 

[D,,D,]{FoG) = {[D,,D,]F)oG + Fo{[Di,Dj]G). (44) 

To show this, use the derivation property of Dj and then the more general result (j43p for Di . 

DiDj{FoG) = D,{DjF oG) + D,{F o DjG) 
[D,DjiFoG)]i = [DijFoG]i + [DiFoDjG]i+ ^ (_i)#»#(/i)+7(/;/i,/2) 

IlUl2=I 

x[DjF]i,[D,G]j,+ (-l)*^')#^'^)+^^'''^''^^i^/i[A,G']/,. (45) 

hUl2=I 

In the third term #(/i) = and in the last term, #(/i) = 0, thus 

Dij{FoG) = D,jF oG + F o DijG + DiF o DjG + {-l)*'^'^*^^^DjF o DiG 
Dji{FoG) = DjiFoG + FoDj,G + DjFoD,G+{-l)*^^*^^^DiFoDjG. (46) 

Combining these we find that the graded commutator of left annihilations is a derivation 

[A, Dj]{F o G) = DiD,{F o G) - (-1)#«#(^')Z),- A(F o G) 

= {D^jF - (-l)#«#(j')Z)j,F) oG + Fo {DijG - (-l)#»#(^^)l),-,G) 

= {[Di,Dj]F)oG + Fo{[D,,Dj\G). (47) 

Iterated commutator of gluonic left annihilation: If all have zero ghost numbers, 

then [[Dj, Dj], Dfc] is a derivation of the shuffle product of two tensors of zero ghost number 
each. This is a consequence of the derivation property of and [Di,Dj\ for shuffle products 
of tensors of zero ghost number. 

[[Di,DjlDk]{FoG) = [Di,Dj]{DkFoG + FoDkG) 

-Dk{[Di, Dj]F o G) - Dk{F o [A, Dj]G) (48) 

Here each of the terms within parentheses is shuffle product of tensors of zero ghost number 
since i,j,k,F,G are. Applying the derivation property of Dk and [Di,Dj] again, four of the 
terms cancel out and we get the desired result 

[[A, Dj],Dk]{F o G) = [[A, Dj],D^]F oG + Fo [[A, DjlDk\G. (49) 

8 <S* for Yang-Mills matrix model with ghosts is a derivation of 
shuffle algebra 

Let us now consider the Yang-Mills matrix model with ghosts introduced in ([2]) 

S = \iT G'^AiAj + tr C'-'''Ai\A,,Ak\ - ^ tr {M, A,\{A^, Ai]g'''g^' (50) 
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This differs from the model considered in (|25p by the addition of the quartic term, which however 
involves only gluons. The latter was studied in [12]. Using our results from [12] and section [5] 
we get the loop equations S^G{S,) = G{^)^^G{^) where 5* is 

= C^^Dj + 3C'^'[D,, D,] - ^g'^g^\D,[Du, Di]]. (51) 

5' is a linear combination of Dj , its graded commutators and gluonic iterated commutator 
[Dj[Dj:, Di]] . Each of these was shown to be a derivation of the graded shuffle product in section 
[71 We conclude that the Schwinger-Dyson operator 5* for Yang- Mills matrix models with ghosts 
is a derivation of the graded shuffle product of ghost number zero correlation tensors. Based on 
the bosonic case in [12], one might have suspected that Di and 5* would, at best, be graded 
derivations upon including ghosts (i.e. satisfy the Leibnitz rule up to a sign). But they turn out 
to be ordinary derivations of the graded shuffle product, since correlations with non-zero ghost 
number vanish in the large- limit. 

9 Discussion 

Physically, the Schwinger-Dyson operator 5* arises from the variation of the action and is 
therefore a classical {h = 0) concept. Indeed, the correlations in the limit h ^ are annihilated 
by the Schwinger-Dyson operator. For Yang-Mills matrix models with ghosts, <S* = C^^Dj + 
3C^^''[Dk, Dj] — ^g^'^g^''[Dj[Dk, Di]] is a linear combination of iterated (anti) commutators of 
the left annihilation operator. The quantum effects in the large- limit are encoded in the 
variation of the measure, or the quadratic term in the loop equations 5*G(^) = G(^)^*G(^) . 
^* are external sources and G(^) is the generating series of correlations, the product on the 
rhs is concatenation. i,j,k label the matrices and in the continuum limit would label space- 
time points. We have shown that there is a mismatch between the 'classical' Schwinger-Dyson 
operator and the 'quantum' concatenation product, so that 5* is not a differential operator in 
the loop equations. This 'mismatch' makes the equations both interesting and hard to solve. 
However, we identified a new commutative shuffle product, with respect to which 5* is indeed 
a differential operator. Shuffle product is the point-wise product of Wilson loops, written in 
terms of correlations. As suggested in jl2] . we can expand concatenation around shuffle to 
exploit the derivation property of left annihilation Di. At O*'^ order, the loop equations become 
differential equations in the graded-shuffle algebra iS*G(^) = G(^) o ^* o G(^) . Moreover, we 
showed that 5* is not just a differential operator with respect to shuffle, but behaves as a first 
order differential operator. This derivation property of 5* leads to a further simplification. First 
we define the reciprocal of G(^) with respect to the graded shuffle product. Since Go = 1 is 
non- vanishing, G(^) has a right reciprocal G(^) o F(^) = 1. Moreover, non-vanishing G/ have 
zero ghost number, so the same is true of the Fj . Thus, the graded shuffle product G o F 
is commutative and the reciprocal is unique and two-sided. Explicitly, i^o = li Fi = —Gi, 
Fij = -Gij + GiGj{l + (-l)#(0#(i)} etc. More generally, for |/| > 1, 

Fi = -Gj- {-iy'^'''''''^Fi,Gi, (52) 

IlUl2=I 

expresses Fj in terms of Gk 's and lower order Fj 's. Iterating, we can find the shuffle reciprocal. 
Now, since 5* satisfies the Leibnitz rule and 5*(1) = 0, 

= S'{G{OoF{C)) =S'GoF + GoS'F S'G = -G o S'F o G. (53) 
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Thus the loop equations become hnear equations for the graded shuffle reciprocal 

S'F = -e. (54) 

This substantial simplification is absent for a generic matrix model whose Schwinger-Dyson op- 
erator is not a derivation of the graded shuffle product. This underscores the potential practical 
importance of the derivation property of the Schwinger-Dyson operator of Yang-Mills matrix 
models, in a possible approximation method based on expanding concatenation around shuffle. 

We can give another physical interpretation. The large- limit is a 'classical' limit since 
U{N) invariants stop fluctuating in this limit, though ?i = 1 is held fixed. However, even this 
large- 'classical' limit is difflcult to solve, partly because we have non-commutative concate- 
nation products left over. Replacing concatenation by commutative shuffle products may be 
thought of as taking a further classical limit. Indeed, when this is done, we found that the 
equations become linear! However, much work still needs to be done. In particular, for some 
models, we found that the loop equations are under-determined [12j. In those cases, the above 
equations have to be supplemented by non-anomalous Ward identities arising from Schwinger- 
Dyson equations that are naively suppressed [18]. We hope to report on more explicit 
calculations in a future publication. 

Mathematically, together with [12], our work shows that it may be fruitful to think of the 
loop equations as living in the differential bi-algebra formed by concatenation, shuffle and their 
derivations. Roughly, the space of based oriented loops (modulo backtracking) on space-time, 
Loop{M) , is to be regarded as a free group on a continuously infinite number of generators 
labeled by the loops 7. The concatenation of loops and reversal of loops are the product and 
inverse operations in this gigantic free group. A typical function on Loop{M) is a Wilson 
loop expectation value with respect to a U{N) connection. The space of functions on this 
Loop{M) is then automatically a commutative Hopf algebra under the commutative point-wise 
product of functions {FG){'y) = F{'y)G{'y) and comultiplication given by concatenation of loops: 
(A'G)(7i,72) = G(7i72). Formally, we also have a dual co-commutative Hopf algebra defined 
via the group algebra of the free group of loops. The Makeenko-Migdal loop equations of large- 
N Yang-Mills theory are defined in this space, using the path derivative of the area derivative 
(the analogue of our Schwinger-Dyson operator), which is a derivation of the pointwise product. 

Hermitian multi-matrix models, regarded as discrete toy-models for gauge theory in physics, 
also provide a finitely generated toy-model for the above mathematical theory on loop space. 
The shuffle product plays the role of the commutative product of functions on Loop{M) . It is 
obtained by expanding the Wilson loop in iterated integrals of gluon correlations. The coproduct 
in the discrete model involves concatenation of tensors (A'(^^) = Jj^^"^®^^) rather than loops. 
One difference (since we are dealing with hermitian rather than unitary matrices), is that we 
do not have inverses for the generators ^* to play the role of reversal of loop orientation. As a 
consequence, the group algebra of the free group of loops is now replaced by the concatenation 
algebra of tensors (free associative algebra), which is the monoid algebra of the free monoid. 
The free associative algebra could also have been obtained as the universal envelope of the free 
Lie algebra. Moreover, left annihilation and its iterated commutators (including the physically 
relevant Schwinger-Dyson operator) form a free Lie algebra of derivations of the shuffle algebra. 
This suggests we should think of the Schwinger-Dyson operator of Yang-Mills theory as a vector 
field on the free group of loops on space time. Thus, despite being a discrete model, our setup 
preserves many of the algebraic and differential structures of the original theory on Loop{M) , 
which is difflcult to study directly. We hope to investigate this correspondence in greater detail. 
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